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1.  INTRODUCTION 


Spike-nosed  configurations  arc  used  for  projectile  applications  against  armored  targets  where 
the  spike  is  used  as  a  stand-off  distance  causing  microseconds  of  advance  time  between  the  time  that 
the  tip  of  the  spike  touches  the  armor  and  the  time  that  the  warhead  (usually  a  shaped  charge) 
detonates.  Spike  nose  configurations  are  also  used  for  a  different  purpose,  namely  reducing  the  drag 
for  blunt  reentry  vehicles  at  hypersonic  speeds  when  drag  and  heating  are  of  major  concern. 

After  World  War  II,  a  new  generation  of  spike-nosed  high  explosive,  anti-tank  (HEAT) 
projectiles  was  developed  in  the  U.S.  and  abroad.  In  the  very  early  stage  of  development,  during  the 
late  1940s  and  early  1950s,  spinning  HEAT  projectiles  were  examined;  but  it  was  quickly  found  that 
spin  reduces  the  depth  of  penetration  in  the  armor.  Therefore,  most  spike-nosed  projectiles,  in  the 
late  1950s  and  thereafter,  were  fm-stabilized  and  were  provided  with  a  tail  boom  and  fins  behind  the 
shaped-charge  warhead. 

Many  experimental  studies  were  made  for  finned,  spike-nosed  projectiles.  Some  wind  tunnel 
tests  can  be  found,  for  example,  in  References  1  through  5.  Some  firing-range  tests  are  reported  in 
References  6  through  9.  The  U.S.  Air  Force,  during  the  1960-70s,  extensively  studied  the  unsteady 
front-shock  flow  phenomenon  (the  buzzing)  for  spike-nosed  reentry  vehicles  at  high  speeds.10-12  Also, 
the  Air  Force  successfully  computed  the  unsteady  buzzing  flow  for  these  reentry  configurations  of 
interest.1315 

At  present,  for  sharp-edged,  spike-nosed  projectiles  of  interest  to  the  Army,  there  is  no  analytic 
or  formal  computational  procedure  that  can  be  systematically  used  to  predict  the  aerodynamics  of  such 
configurations.  The  Army  has  relied  so  far  on  direct  wind  tunnel  tests  followed  by  live  firing  of 
projectiles  in  the  firing  ranges.  The  present  study  was  made  to  establish  such  a  systematic,  numerical, 
predictive  technique.  Therefore,  validation  of  the  predictive  technique  against  range  or  wind  tunnel  data 
is  of  vital  importance  for  assessing  the  numerical  capability.  Although  the  final  objective  is  the 
application  to  Finned,  spikc-noscd  projectiles,  this  study,  being  a  first  step  toward  that  goal,  limited 
itself  to  unflnned,  spike-nosed  configurations  in  an  attempt  to  focus  on  the  spike-nosed  flow  with  its 
complex  features.  These  features  include  dual  flow  modes,  large  separation  regions,  and  unsteadiness. 
By  establishing  that  such  flows  can  be  systematically  and  successfully  computed,  the  doors  will  be 
opened  for  future  work  to  tack'c  similar  configurations  with  added  booms  and  fins. 


The  advances  in  the  zonal  gridding  and  overlapping  techniques  made  this  study  possible  for 
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the  present  sharp-edged  configurations.  This  represents  the  first  known  application  of  computational 
fluid  dynamics  to  Army  spike-nosed  projectile  shapes. 


2.  ABOUT  THE  TEST  CASES 

Very  few  experimental  results  arc  available  for  unfinned,  spike-nosed  projectiles  in  contrast  to 
those  of  finned  ones  for  the  reason  stated  earlier.  For  the  sharp-edged  configuration,  only  the  wind 
tunnel  tests  of  Platou16  arc  applicable.  Also,  some  firing  range  data  arc  available,  but  for  configurations 
with  a  tripping  ring.1'  Very  recent  wind  tunnel  tests  were  made  by  Koenig,  cl  al..18  while  this  work 
was  being  completed,  for  very  similar  configurations  w-ith  different  spike  lengths  at  Mach  nuirl>ers 
between  0.8  to  1.5  and  zero  angle  of  incidence.  Results  of  Reference  18  could  be  used  for  further 
validations  in  the  future. 

Three  cases  from  Reference  16  were  chosen.  They  are  shown  in  Figures  1-3.  The  wind  tunnel 
tests  were  made  in  1950  at  the  Ballistic  Research  Laboratory  supersonic  wind  tunnels  at  M=1.72  for 
angles  -10°  <  a  <  +10°.  The  Reynolds  number  was  4.86x10*  per  foot,  To=10O  F°,  and  P0=1.26 
atmospheres.  The  model  diameter  is  2.5  inch  for  all  three  configurations.  The  first  two  configurations 
were  reported  not  to  have  dual  flow  modes  at  this  Mach  number  and  range  of  a.  The  flow  features 
of  both  modes  are  depicted  in  Figure  4.  The  third  configuration  was  reported  to  have  had  the  dual 
flow  modes  (the  high-drag,  open-flow  mode  and  the  low-drag,  closed  flow  mode).  However,  the  low- 
drag  mode  briefly  occurred  while  increasing  a  and  was  captured  on  a  schliercn  photograph,  but  the 
drag  force  itself  was  not  measured.  The  high-drag  mode  then  persisted,  while  the  lower  drag  mode 
could  never  be  recovered  again  during  the  tests.16 


3.  GOVERNING  EQUATIONS 


The  compressible,  turbulent  Navicr-Stokes  equations  for  axisymmeiric  and  two-dimensional  flow 
can  be  expressed19  in  the  following  strong  conservation  form,  in  which  the  dependent  variables  p,  u, 
v,  and  c  are  mass  averaged,  with  c  being  the  specific  total  energy,  T  being  the  temperature,  p  and  p 
being  mean  density  and  pressure,  respectively,  and  l  denoting  time: 
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where  p.  is  molecular  viscosity,  e  is  the  turbulent  eddy  viscosity,  and  P  =  1  or  0  for  axisymmctric 
and  iwo-dimcnsional  cases,  respectively. 


The  air  is  assumed  to  be  a  perfect  gas,  satisfying  the  equation  of  state  p  =  pRT,  where  R  is 
the  gas  constant  (1,716  ft2/sec2  -  °R  for  air).  For  the  dependence  l  •  aninar  viscosity  on  temperature, 
Sutherland’s  law  was  used: 


p  =  2.270 


T  +  198.6 


10 8 
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The  laminar  and  turbulent  Prandtl  numbers,  Pr  and  Pr„  were  assumed  constant  with  values  of 
0.72  and  0.9,  respectively.  The  ratio  of  specific  heats,  y,  was  also  assumed  constant  and  equal  to  1.4. 
C,  and  Cp  arc  specific  heal  capacities  at  constant  volume  and  constant  pressure,  respectively: 

Cv  =  4290  ftVsec2  -  °R. 
and 

Cp  =  6006  ft2/scc2  -  °R  for  air. 


The  total  energy  per  unit  mass,  e,  is  given  by: 

c  =  C,T  +  (1/2)  (u2+  v2). 


In  the  5,  -  r\  computational  plane,  Equations  la  and  b  are  transformed  to  the  conservation  law 
form,  and  the  equations  can  he  found,  for  example,  in  Reference  19. 

31  Turbulence  Model.  Turbulence  is  modeled  through  the  algebraic  eddy-viscosity  model 
of  Baldwin  and  Lomax.20  This  model  employs  the  two-layer  concept  (inner  and  outer).  The  inner 
layer  is  near  the  walls  and  is  modeled  as: 
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The  distance  normal  to  the  surface  is  y;  A"  =  26;  k  =  0.40  is  the  von  Karman  constant;  and  the 
subscript  w  denotes  values  at  the  surface. 

The  model  switches  from  the  inner  to  the  outer  region  at  the  smallest  value  of  y  for  which  the 
inner  and  outer  values  of  the  eddy-viscosity  arc  equal  (i.c.,  e,  =  e0).  The  e  for  the  outer  layer  is  given 
by 
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The  value  of  y  at  which  FinM  occurs  is  yra„. 
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Due  to  the  perpendicular  surfaces  of  the  spike  surfaces  at  the  nose  tip  and  at  the  facing  shouldu,  the 
normal  distance  to  the  wall,  y,  in  Equation  3  is  difficult  to  assign.21  This  problem  was  solved  in 
Reference  21  by  measuring  the  y  along  a  45°  ray  emanating  from  the  point  of  intersection  of  the  two 
perpendicular  walls. 


4.  ABOUT  THE  CODE,  GRID,  AND  COMPUTATIONS 

4.1  The  Code.  The  code  was  developed  by  Patel  and  Sturck19.  It  utilizes  the  familiar  and 
robust,  explicit,  timc-dcpcndcnt  method  of  McCormack.  The  code  was  vectorized  and  is  run  on  a 
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Cray-XMF/48  Machine.  The  present  computations  were  all  run  in  serial  arithmetic  mode.  The  zonal 
grid  and  overlap  provided  in  the  code  are  represented  by  eight  different  available  zones  (which  can  be 
increased  if  so  desired).  The  user  prescribes  the  overlap  between  regions  along  one  line  of  adjacent 
zones  (interface).  A  global,  uniform,  time  step  was  used  herein  against  grid-varying  time  steps  to 
simulate  time-accurate  solutions.  The  time  step  is  determined  from  the  Courant-Fredric.k-Levy  (CFL) 
condition,  with  a  factor  of  about  0.6  being  used  as  the  Courant  number. 

4.2  Boundary  Conditions.  No-slip  conditions  are  specified  on  all  wall  surfaces.  The  incoming 
flow  conditions  are  assumed  to  be  of  uniform  profiles  with  free-strcam  temperature  of  T_  =  520  R°, 
p_  =  14.7  psi,  =  1.72. 

The  outgoing  conditions  at  the  end  of  the  projectile  were  imposed  as  zero  gradients  parallel  to 
the  body  axis  direction. 

The  outer  boundary  conditions  were  imposed  as  nonrcflective  conditions,  i.c.,  zero-gradient 
conditions  along  characteristic  lines  for  all  variables.  The  characteristic  direction  is  determined  from 
the  local  velocity  and  temperature.  This  approach  allows  setting  the  "outer"  field  close  to  the  body 
without  the  penalty  of  any  unnecessary  approximations  regarding  shock  reflection,  or  zero-gradient 
conditions. 

At  the  symmetry  line,  ahead  of  the  spike  tip,  a  two-point,  zero-gradient,  boundary  condition 
is  imposed  on  the  solved  variables. 

4.3  Initial  Conditions.  Computations  were  started  using  free-strcam  values  everywhere  in  the 
domain.  These  values  are  for  free  stream  velocity,  pressure,  and  temperature.  The  density  and  specific 
total  energy  are  computed  accordingly,  using  the  equation  of  state  and  the  definition  of  the  specific 
total  energy. 

44  The  Grid.  Three  different  grid  zones  were  used  in  the  computation.  Those  zones  and  the 
extent  of  the  computational  domain  arc  depicted  in  Figure  5. 

For  the  first  configuration,  the  grids  used  for  the  three  zones  arc  (15x48),  (26x39),  and  (21x25), 
respectively.  The  first  and  second  arguments  in  the  parentheses  refer  to  the  axial  and  radial  directions, 
respectively.  This  grid  has  2,259  total  points  and  is  equivalent  to  a  (48x48)  grid. 
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One  restriction  in  the  present  grid  overlapping  technique  is  the  requirement  that  no  interpolation 
is  allowed  at  the  interface  line  between  zones.  Thus,  each  point  on  either  side  of  any  two  zones  must 
have  exactly  the  same  coordinates.  This  restriction  represents  some  constraint  in  the  flexibility  of  the 
grid  distribution  and  may  be  alleviated  in  future  development  of  the  code.  Meanwhile,  to  accommodate 
this  restriction,  one  has  to  accept  unnecessary  clustering  of  points  in  some  locations.  Figure  6  shows 
the  clustered  points  along  lines  parallel  to  the  top  body  surface,  where  clustering  is  needed  near  the 
body  to  resolve  the  turbulent  boundary  layer.  Figure  6  shows  the  overall  grid  distribution  for  the  first 
projectile  configuration. 

The  sharp  cone  spike  configuration  was  also  modeled  using  three  zones,  but  with  grid  sizes 
of  (15x39),  (26x39).  and  (21x25).  respectively.  This  grid  totals  2,124  points,  or  the  equivalent  of  a 
(46x46  mesh),  and  is  depicted  in  Figure  7. 

The  third  configuration  was  computed  using  zones  of  sizes  (15x44),  (21x35),  and  (31x35). 
This  grid  totals  2,480  points  or  (50x50),  approximately.  This  grid  for  configuration  3  is  given  in 
Figure  8.  The  grids  were  generated  using  a  simple  algebraic  (exponential)  formula  in  each  zone  in 
both  the  axial  and  radial  directions. 


5.  RESULTS 

5.1  Configuration  1.  The  low-drag  mode  was  obtained  by  straightforward  computation, 
assuming  the  flow  to  be  fully  turbulent  everywhere.  Unsteadiness  occurred  in  the  computation,  but, 
when  the  turbulence  level  was  reduced  to  0.2  of  its  value  at  each  point,  the  flow  became  steady. 
The  Mach  contours  arc  provided  in  Figure  9.  Comparison  with  the  schlicrcn  photograph  of  Reference 
16  indicated  good  agreement  of  flow  features.  The  bow  shock  stand-off  distance  was  0.5  d5  as 
theoretically  predicted,  and  the  bow  shock  angle  away  from  the  body  nose  tip  was  aboui  40°,  the  same 
as  can  be  determined  from  the  schlicrcn.  The  computation  converged  satisfactorily  after  4,000  time 
steps,  although  the  code  was  later  run  to  12,000  steps  to  assure  the  stability  of  computations.  The 
computer  CPU  time  was  40  minutes  on  the  Cray-XMP/48  for  the  4,000  steps. 

Several  numerical  experiments  were  made  to  obtain  the  high-drag  mode,  even  though  that 
mode  was  not  confirmed  during  the  tests  of  Reference  16.  The  high-drag  mode  was  easily  obtained 
by  freezing  (i.c.,  not  updating)  the  turbulence  level  after  1,000  steps  and  by  simultaneously  relaxing 
the  time  step  at  each  point  to  0.6  of  its  local  Courant  time  step  value.  The  solution  also  converged 
satisfactorily  after  6.000  time  steps,  using  60  minutes  on  the  same  computer.  Figure  10  depicts  the 
flow  field  as  presented  by  the  Mach  contours  for  this  high-drag  mode. 
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Comparison  between  the  local  flow  fields  of  the  two  modes  near  the  spike  tip  is  given  in 
Figures  11  and  12,  respectively.  For  the  high-drag  mode,  there  is  an  expansion  fan  near  the  tip,  which 
is  followed  immediately  by  coalescence  of  compression  waves  facing  the  separation  region.  These 
compression  waves  coalesce  into  a  shock  which  faces  the  facing  shoulder  of  the  projectile,  thus  raising 
the  pressure  behind  it  and  also  that  on  the  facing  wall.  'Hus  higher  pressure  results  in  the  higher  drag 
of  the  projectile.  Surprisingly,  the  corresponding  flow  detail  near  the  facing  shoulder  differs  very 
slightly  for  the  two  modes.  It  was  expected  that  larger  differences  would  be  observed  there.  Figures 
13  and  14  provide  the  details  for  those  modes  near  the  facing  shoulder. 

The  forcbody  drag  coefficient  for  the  computed  geometry  (low-drag)  was  0.337.  The  drag  due 
to  the  rotating  band,  which  is  shown  in  Figures  1-3  but  was  not  modeled  in  the  computation,  was 
estimated22  to  be  2%  of  the  total  drag  at  M  =  1.72.  The  computed  drag  is,  therefore,  provided  as 
0  344,  while  the  wind  tunnel  measurement  given  in  Reference  16  is  0.351  for  the  forcbody  drag. 
Reference  16  provides  the  nei  forcbody  drag  without  any  reference  to  base  drag  corrections.  The 
computation,  thus,  undcrprcdicts  the  measurement  by  2%.  Considering  the  tunnel  measurement 
accuracy,  one  can  conclude  that  these  flow  results  are  very  assuring  and  useful.  This  result  is  shown 
in  Figure  15.  For  the  high-drag  mode,  the  computed  drag  coefficient  was  0.402  and  is  provided  as 
0.410  when  including  the  2%  rotating  band  effect.  The  high-drag  mode,  therefore,  resulted  in  a 
19%  increase  in  drag  over  the  low-drag  mode. 

5.2  Configuration  2.  The  high-drag  mode  was  obtained  first  when  the  computation  of 
Configuration  1  was  repeated,  assuming  fully  turbulent  flow  everywhere  and  no  reduction  in  the 
turbulence  level.  The  computation  was  slower  in  converging,  requiring  7,000  time  steps  for  satisfactory 
convergence. 

To  obtain  the  low-drag  mode,  which  is  the  mode  reported16  to  occur,  laminar  flow  and 
transition  should  be  allowed  to  occur  on  the  cone.  Therefore,  laminar  flow  was  allowed  on  the  cone, 
and  transition  was  allowed  to  occur  only  along  the  middle  third  of  the  whole  spike  length.  This  was 
based  on  estimates  of  location  of  transition  Oocal  Reynolds  number),  which  were  evaluated  using 
References  23  and  24.  The  low-drag  mode  was  immediately  obtained,  but  with  slower  convergence 
rate.  Convergence  required  about  12,000  time  steps,  requiring  120  minutes  on  the  Cray-XMP/48 
Machine. 

The  two  How  modes  are  depicted  in  Figures  16  and  17.  In  Figure  16,  the  flow  seems  to  slide 
over  the  separated  region  of  the  spike,  while  in  Figure  17,  for  the  high-drag  mode,  there  is  a 
compression  wave  appearing  at  the  beginning  of  the  separated  region,  thus  signifying  flow  path  turning. 
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The  forebody  drag  coefficient  for  the  low-drag  case  was  0.314  (including  the  2%  rotating 
band  drag),  compared  to  0.321  for  the  wind  tunnel  measurement.  Again,  computation  is  within  2% 
of  the  measurement.  Figure  18  depicts  the  comparison  in  the  drag  value.  The  high-drag  mode  drag 
coefficient  was  computed  as  0.397  (including  the  2%)  and,  thus,  is  26%  higher  than  that  of  the  low- 
drag  mode.  The  high  drag  mode  was  not  observed  during  the  tests  of  Reference  16. 

5.3  Confieuration  3.  The  high-drag  mode  was  obtained  readily  when  the  computation,  similar 
to  that  for  Configuration  1,  was  applied  here.  Fully  turbulent  flow  everywhere  with  no  reduction  in 
turbulence  level  was  applied.  Computations  required  only  4,000  steps  for  satisfactory  convergence. 

Several  numerical  attempts  were  made  to  obtain  the  low-drag  mode.  It  was  found  that,  by 
imposing  a  3.5%  cross  flow  (i.e.,  v  =  0.035  VJ  in  the  free  stream  in  zone  1,  the  flow  mode  was 
readily  obtained.  This  3.5%  cross  flow  falsely  simulated  a  pseudo  2°  angle  of  attack.  Although  this 
is  not  truly  an  angle-of-attack  effect,  the  cross  velocity  is  an  influence  that  can  be  related  to  an  angle 
of  attack. 

These  flow  fields  arc  depicted  in  Figures  19  and  20  for  low-drag  and  high-urag  modes, 
respectively.  An  excellent  agreement  for  the  high-drag  mode  was  obtained  with  the  schlieren 
photograph  of  Reference  16.  One  feature  is  the  existence  of  a  "kink"  in  the  compression  wave,  which 
emanates  from  the  impact  of  the  flow  with  the  separation  region.  This  kink  was  questioned  at  first, 
but,  when  the  schlieren  photo  had  been  examined  carefully,  the  kink  was  found  easily.  Also,  all  the 
shock  and  expansion  wave  angles  (away  from  the  body)  were  found  to  be  within  4°  of  the  values 
measured  from  the  schlieren  photograph  of  Reference  16. 

The  forebody  drag  coefficient  for  the  high-drag  mode  was  computed  as  0.478  (including  the 
2%  increase  due  to  the  rotating  band)  compared  to  0.306  for  the  low-drag  case.  This  represents  a 
56%  increase  in  drag.  The  wind  tunnel  measurement  for  the  high-drag  mode  was  reported'6  to  be 
0.555.  This  large  discrepancy  between  the  computed  and  measured  values  is  still  unresolved. 
However,  when  examining  this  particular  case  in  Reference  16,  the  value  of  0.555  seems  to  be 
particularly  high  in  comparison  with  the  remaining  cases  tested.  The  data  of  Reference  16  for  this  case 
were  presented  with  only  one  point  on  one  figure  in  the  report.  There  arc  no  cross-checkcd  values  or 
any  tabulated  results  for  positive  verification  of  this  value.  Table  1  provides  a  summary'  of  all  obtained 
results  and  a  comparison  with  test  data. 

5.4  Gnd  Size  Effect.  To  examine  the  result  obtained  for  the  high-drag  mode  i  .figuration 
3,  the  computations  were  performed  again  using  a  larger  number  of  points  to  assess  whether  that 
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discrepancy  was  due  to  inadequate  grid  size.  Configuration  3  was  computed  first  using  zones  with 
(15x44).  (21x35).  and  (31x35)  points.  This  grid  totals  2,480  points,  with  a  (50x50)  mesh  equivalence. 
The  grid  was  then  increased  to  (15x54),  (31x45),  and  (41x45),  thus  totalling  4,050  points,  which  is 
equivalent  to  a  (63x63)  mesh.  The  drag  coefficient  changed  from  0.4621  to  0.4690,  a  change  of  only 
1.5%.  Therefore,  it  was  assumed  that  the  grid  size  is  appropriate  for  most  purposes. 


6.  SUMMARY  AND  CONCLUSIONS 

Three  different  spike-nosed  projectile  configurations  were  computed  at  Mach  =  1.72,  and  the 
results  are  compared  to  wind  tunnel  measurements.  The  computed  drag  coefficients  are  in  very  good 
agreement  with  the  measured  values.  Computed  values  arc  within  2%  of  the  measurements,  which  is 
within  the  diag-measurement  accuracy  itself.  The  high-drag  mode  computed  for  the  third  configuration 
provided  a  considerably  lower  drag  than  the  measured  value,  although  the  detailed  flow  features 
compared  rather  accurately  with  the  schlicrcn  photograph  of  the  test.  Because  confidence  was  gained 
from  the  two  previous  computed  cases,  it  is  believed  that  the  drag  measurement  for  this  particular  case 
is  quite  high  and  may  also  be  in  error.  This  belief  is  supported  by  observing  the  results  of  20  similar 
spike-nosed  configurations  tested  during  the  same  test  period.  These  drag  data  are  also  reported  in 
Reference  16. 

Two  interesting  obstacles  are  faced  in  this  study  which  cannot  be  resolved  decisively.  First, 
the  possibility  of  two  flow  modes  always  exists.  'I'herefore,  one  cannot  determine,  a  priori,  which 
mode  the  numerical  procedure  will  favor.  Also,  it  is  not  known,  a  priori,  which  one  will  physically 
occur  at  particular  wind  tunnel  conditions  or  under  free-flighl,  firing-range  conditions.  The  second 
obstacle  faced  was  the  strong  role  of  the  turbulence  eddy  viscosity  value  level  and  model  on 
influencing  the  computation  (possibly  due  to  the  large  separation  region)  toward  one  particular  flow 
mode.  It  was  not  determined  whether  this  influence  is  purely  a  numerical  problem  or  it  has  a  parallel 
in  nature  where  flow  turbulence  in  the  tunnel  or  in  free  flight  may  trigger  a  particular  flow  mode. 

Finally,  this  study  provided  a  straightforward  and  systematic  capability  for  computing  such 
difficult  configurations.  The  present  work  represents  an  advance  in  the  application  of  CFD  techniques. 
The  computations,  in  addition,  have  provided  dual  flow  modes  where  the  wind  tunnel  experiment  had 
only  revealed  one  mode  under  certain  tunnel  conditions.  It  is  not  known  if  some  of  these  computed 
dual  modes  arc  superficial  or  whether  real-life  tests  had  favored  only  one  mode  which  is  more 
dominant.  Thus,  these  computations  may  spur  the  need  for  extensive  and  delicate  variations  in  test 
conditions  to  verify  the  existence  of  these  modes  at  these  flow  conditions. 
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Wind  tunnel  model  of  tne  second  projectile  confi 


Wind  tunnel  model  of  the  third  projectile  confitnjration. 


Figure  4.  Typical  flow  features  for  the  low-  and  hi 
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Figure  7.  Zonal  grid  for  ihe  second  projectile  configuration. 


GRID 


’8 


Figure  8.  Zonal  grid  for  ihe  third  pnojeclile  configuration. 
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Figure  9.  Mach  number  contours  for  the  low-drag  mode  for  Con fipu ration  1 
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Figure  11.  Dclails  of  Mach  number  contou 
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Figure  12.  Details  of  Mach  number  contours  near  the  spike  tin  (higMrag.  Configuration  1 
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Figure  14.  Deiails  of  Mach  number  contours  near  ihe  facing  shoulder  (hiph-i 
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Figure  15.  Forcbodv  drag  coefficient  comparison  for  Configuration 
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Mach  number  contours  for  Configuration  2  (low-dra; 
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Figure  17.  Mach  number  contours  for  Configure! 
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Figure  18.  Fprebody  drag  coefficient  comparison  for  Configuration  2  flow-drag  modeV 
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Figure  19.  Mach  number  contours  for  Configuration  3  down 
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TABLE  1.  Comparison  of  Results  at  Mach  =  1.72  . 


Configuration 

C, 

Low-Drae  Mode 

Hieh-Drae  Mode 

Wind 

Tunnel 

Present 

Computations* 

Wind 

Tunnel 

Present 

Computations* 

1 

0,351 

0.344 

_b 

0.410 

2 

0.321 

0.314 

b 

0.395 

3 

C 

0.306 

0.555 

0.469 

“These  values  include  an  added  2%  due  to  rotating  band  pressure  drag,22  at 
M  -  1.72  . 


‘This  flow  mode  was  not  reported  in  the  wind  tunnel  experiment.16 

This  mode  was  observed  and  reported  only  in  a  schlicrcn  photograph  but 
quickly  disappeared  and  could  not  be  recovered16  in  the  wind  tunnel  for  actual 
measurement. 
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LIST  OF  SYMBOLS 


Nomenclature 


case 


Ajtf  =  reference  area,  (7td2/4) 

CD  =  drag  coefficient,  drag  force/(.5  p_  V2_  A,d) 

Cp  =  specific  heat  under  constant  pressure 
Cv  =  specific  heat  under  constant  volume 
d  =  reference  diameter 
d,  s  spike  diameter 
e  =  specific  total  energy 
M  =  Mach  number 
p  =  static  pressure 
Re  =  Reynolds  number 

u.v  =  velocity  components  in  the  x.y  directions 
V„  =  free  stream  velocity 

x,y  =  Cancsian  coordinates  for  2-D  case,  axial  and  radial  coordinates  for  axisymmctric 


Creek  symbols 

a  =  angle  of  attack 

7  =  ratio  of  specific  heats  for  air 

p  =  density 

p  =  laminar  (molecular)  viscosity  coefficient 

£  =  turbulent  eddy  viscosity  coefficient 

4.11  =  transformed  coordinates  in  the  computational  plane  for  the  coordinates  x,y 


Subscripts 

o  =  denotes  stagnation  (total)  condition 
00  =  free  stream  condition 
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